Feasibility of a storage ring for polar molecules in strong-field-seeking states 
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We show, through modehng and simulation, that it is feasible to construct a storage ring that will 
store dense bunches of strong-field-seeking polar molecules at 30 m/s (kinetic energy of 2K) and hold 
them, for several minutes, against losses due to defocusing, oscillations, and diffusion. The ring, 3 m 
in diameter, has straight sections that afford access to the stored molecules and a lattice structure 
that may be adapted for evaporative cooling. Simulation is done using a newly-developed code that 
tracks the particles, in time, through 400 turns; it accounts for longitudinal velocity changes as a 
function of external electric field, focusing and deflection nonlinearities, and the effects of gravity. 
An injector, decelerator, and source are included and intensities are calculated. 

PACS numbers: 29.20 Dh, 41.75.Lx, 33.80.Ps, 39.90.-fd, 33. 55. Be 



I. INTRODUCTION 

To date, all gaseous quantum condensates have been 
produced by evaporative cooling of confined atoms. Con- 
finement is necessary to thermally isolate the parti- 
cles from the warmer environment and long confinement 
times are necessary because the evaporative cooling pro- 
cess can take tens of seconds. 

Strong magnetic field gradients have been used to 
confine neutral paramagnetic molecules jlj and electric- 
field gradients have been used to confine neutral polar 
molecules in electrostatic traps and in toroidal stor- 
age rings "s", In addition, polar molecule confinement 
in a synchrotron storage ring has been modeled 0. 

All of these methods use molecules or atoms in weak- 
field-seeking states, whose binding energy decreases in 
the field. These states are not the lowest energy state and 
are therefore subject to coUisional relaxation. In alkali 
atoms, the relaxation rates from the stretched hyperfine 
levels (mp = F) is small. But in magnetically trapped 
paramagnetic molecules Jyi and in electrically confined 
polar molecules 0, S IS ^| , the relaxation rate can be 
large enough to prevent achieving the confinement time 
needed for evaporative cooling. 

CoUisional relaxation will be absent for polar molecules 
in their lowest rotational state. This ground state is 
strong-field-seeking, as are all rotational states in the 
limit of strong electric field. The technical challenges 
of storing molecules in a strong-field-seeking state have 
not been previously addressed. The major challenge is 
focusing these molecules because electrostatic lenses can 



focus strong-field-seeking molecules in only one trans- 
verse plane while defocusing in the other. Therefore 
alternating-gradient focusing is required. 

For experiments on molecules in strong-field-seeking 
states, a storage ring has some useful features not gener- 
ally found in traps. The ring has a beam geometry with 
field- free regions accessible to experiments, and it can si- 
multaneously store many bunches of particles producing 
a large flux of molecules. 

In this paper we show, by modeling and simulation, 
that it is feasible to construct a storage ring (Fig^ that 
will store a symmetric-top molecule (methyl fluoride) in 
the J = state, at a kinetic energy of 2 K (30 m/s), 
and by extension other molecules and velocities. In the 
storage ring, bunching electrodes hold the molecules in 
a string of short bunches. The molecules are calculated 
to be stable against losses due to defocusing, oscillations, 
and diffusion for over two minutes. We also model a 
decelerator for slowing the molecules to 30 m/s, and an 
injector for loading the storage ring. 

A storage ring in which the density of the molecules in 
a bunch is allowed to vary around the ring, can provide a 
mechanism for evaporative cooling. Regions of high den- 
sity speed the thermalization of the molecules. In regions 
of low density the molecules can become spatially sepa- 
rated due to their velocity spread allowing the hottest 
molecules to be removed. 



II. FORCES DUE TO ELECTRIC FIELD 
GRADIENTS 

A. Focusing and Deflection Using Multipole Fields 
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A brief description of focusing and deflecting a beam 
of molecules using electrostatic multipole fields is given 
below. Additional details of beam transport and focusing 
of molecules in strong-field-seeking states, with specific 
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FIG. 1: Layout of the storage ring. Each octant contains a 
buncher and a pair of alternating-gradient focusing triplets to 
match the beam traversing from the straight sections to the 
bend sections. A bend section contains combined bend and al- 
ternating gradient focusing elements. The focusing and bend 
elements have time-independent electric fields. An injection 
line is located in one of the straight sections 



application to methyl fluoride in the J = state, may be 
found in Kalnins et al. (llj . 

The guide field in a storage ring for molecules in strong- 
field-seeking states must provide all the functions, such 
as focusing, bending, and bunching, that are used in a 
ring for charged particles but with forces that arise from 
gradients of the magnitude of the electric field. 

In a pure quadrupole or sextupole field, the total elec- 
tric field increases radially and the force on a molecule, 
in a strong-field-seeking state, is away from the center- 
line in all transverse directions. Therefore a dipole com- 
ponent must be added to remove the double-defocusing, 
and obtain focusing in one transverse direction while still 
defocusing in the other. The force on a molecule is given 
by the gradient of its Stark potential energy, W{E): 



tipole potential used to bend and focus a molecule is 

* = Eo[y + A^xy + Asix^y - ^y^)] 



(2) 



where Eq is the dipole field strength, and A2 and A3 
are the relative quadrupole and sextupole component 
strengths. 

For the Stark energy in the high-field limit, the forces 
to second order are: 



E, 



deEo[A2 

deEoHA 



2Asx~^{Al-AA,As)y'] 
- 2^3)2/ - (^2 - 4A2A3)xy] (3) 



We see that a combined dipole and sextupole (A^) field 
lens will focus in one plane, while defocusing in the other. 
To deflect the molecule we must add a quadrupole {A2) 
component. This also defocuses the beam in the y di- 
rection and stronger sextupole (A3) strengths are needed 

To obtain net focusing in both transverse planes, the 
lenses are arranged in a sequence with gradients alter- 
nating in sign {A3 < for x-focusing and A3 >0 for 
y- focusing). 



B. Other Effects 

When a molecule in a strong-fleld-seeking state enters 
the field of an electrode pair it is accelerated longitudi- 
nally, and upon exiting the field it is decelerated. Also, 
the fringing field is stronger away from the midplane and 
this causes a net defocusing force in the direction of the 
electric field. Between successive sets of electrodes, this 
unwanted defocusing is reduced if the dipole fields are of 
the same polarity and strength. 

Longitudinal bunching, as in a charged-particle ring, 
requires a pulsed field. The field is ramped in a sawtooth 
or sine- wave form and the time-dependent acceleration is 
the net difference between the fields when entering and 
when exiting. 

The effect of gravity is small but not negligible for 30 
m/s molecules in this ring. The vertical orbit will be 
distorted and an orbit correction must be applied. 



C. Equations of Motion 



F = ~\/W{E) = -—VE (1) 
dE 

where E is the magnitude of an external field. 

The Stark energy of the molecular level is in general 
a nonlinear function and is described for methyl fluoride 
in the J = rotational state in Ref. j^J- the limit 
of large E, W{E) —deE where de is the molecule's 
electric dipole moment. 

The transverse {x horizontal, y vertical) electric mul- 



The equations of motion of a molecule in the ring are 
obtained from the Hamiltonian: 

H = Ho + W{E) - gy (4) 

where W{E) is the Stark energy, g is the acceleration due 
to gravity, and Hq is the kinetic energy which in a bend 
region is: 
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where Px and Py are the transverse momenta, Pg is the 
angular momentum and p is the bend radius. In straight 
sections the last term is replaced by the square of the 
longitudinal momentum, Pz^- 

The longitudinal variation of the Stark energy at the 
ends of electrodes (treated here as a step function) adds 
or subtracts from the kinetic energy, the change in lon- 
gitudinal velocity being about ±10%. 

Vertical defocusing in a fringe field is derived from the 
longitudinal variation of the field on the midplane and to 
lowest order is: 

iF)r- --^ , (6) 




Time 



FIG. 2: A molecule at the bunch center enters and exits the 
buncher when the field is the same and receives no net accel- 
eration. For a molecule that arrives later, the entering field 
is stronger than at its exit; it is accelerated and it then drifts 
downstream toward the bunch center. 



III. STORAGE RING DESIGN 
A. Molecule and Energy 

The principles and techniques we use apply to all po- 
lar molecules in strong-field-seeking states. We choose 
methyl fiuoride (CH3F) as our reference molecule because 
it is a nearly symmetric rotor with a large electric dipole 
moment of 4 = 6.2 x 10"^'' C-m (1.84 D). It has a mod- 
erate rotational constant oi B = 0.88 cm~^ and a simple 
level structure with a J = = rotational ground state. 
The rotational constant is large enough to limit the num- 
ber of rotational levels populated in the beam from a jet- 
source but still small enough to allow for a large Stark 
effect at moderate electric fields. Methyl fluoride is also 
a gas at room temperature. 

The velocity of 30 m/s (kinetic energy of about 2K) is 
low enough to make for a compact ring, yet keep small 
the effects of gravity. 

B. Ring Lattice 

Long straight regions free of focusing electrodes make 
the stored beam accessible for experiments and give space 
for injection and extraction. Molecules, in order to drift 
through the straight section without loss, must have only 
small divergences and therefore a large beam width. In a 
bending region, we need strong deflecting forces to min- 
imize the bend radius for overall compactness. These 
strong forces call for a small beam width to avoid non- 
linearities. To make the transition (match) from straight 
sections to arc sections, triplets (Ql, Q2, Q3) of focusing 
lenses are placed at the ends of the straight sections, as 
shown in Fig. ^ 

In each of the eight bend regions, there are five elec- 
trode pairs; each has a combined dipole and quadrupole 
field to provide the strong deflecting force. To this is 
added a sextupole component, the gradient of which al- 
ternates in sign. 

The electrode parameters are given in Table where 
Q are focusing elements and BE and BD are combined 



bend and focusing elements. Each arc is a series of BE 
and BD elements: iBE+BD+BE-|-BD-hiBE. 

In this sequence of lenses with alternating gradients, 
the molecules execute oscillatory transverse motions. 
The parameters of BE and BD are chosen such that the 
phases of these horizontal and vertical motions each ad- 
vance through an angle of 27r in each octant of arc. The 
parameters of Ql, Q2 and Q3 are varied to find values 
that produce large dynamic aperture and momentum ac- 
ceptance. The decapole coefficient, of Q2, which adds 
the term EoA^{x^y-2x'^y'^ + \y^) to the potential Eq©, 
is introduced to reduce the nonlinearity of Q2 focusing 
where the beam is at it's largest. Eor longitudinal con- 
finement with many short bunches, we use eight bunchers 
in the ring; each has a short uniform field that is pulsed 
in time as illustrated in Eig. [3 

Molecules with different energies have their closed or- 
bits radially separated in the arcs and perhaps elsewhere 
in the ring. If this dispersion of orbits is present at a 
buncher, the energy change from the buncher produces 
a shift in the orbit and an increment in the radial os- 
cillation. This is called synchro-betatron coupling and 
to avoid growth of radial oscillation amplitude, the dis- 
persion of orbits must be made zero at the bunchers. 
With the phases of the vertical and horizontal motions 
advancing through an angle of 27r in each octant, as noted 
above, the dispersion becomes zero at all eight buncher 
locations. 



TABLE I: Parameters of Storage Ring Electrodes 





Eo 


L 


A2 


Ai 


As 




(MV/m) 


(cm) 


(m-2) 


(m-^) 


(m-^) 


Ql 


3.0 


3.34 





2000 





Q2 


4.0 


3.71 





-2000 


-1.28x10*^ 


Q3 


4.0 


2.85 





2000 





BF 


7.85 


4.00 


-10.55 


-2296 





BD 


7.85 


4.00 


-10.55 


2343 
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C. Numerical Modeling and Simulation 

The lattice parameters (Table ^) are found by numeri- 
cal calculations using a newly-developed simulation code 
that tracks the particles in time (rather than in longi- 
tudinal position) to account for the longitudinal velocity 
changes as a function of the external field. The tracking 
code includes the effects of nonlinearities, gravity and 
the longitudinal kick at the bunchers. The effect of each 
fringe field (Eq. (S)) in every element has been integrated 
and replaced by a vertically defocusing thin lens. The pa- 
rameters in Table result in the ring performance listed 
in Table. ITU and shown in Figures 13 and 0| 

TABLE II: Ring Parameters 



Parameter 


Value 


Circumference (m) 


9.850 


Circulation period (s) 


0.3121 


Velocity in free space (m/s) 


30.0 


Symmetry of the ring 


8 


Bending radius (m) 


0.60 


Long straight section (m) 


0.40 


Beta function* (m) 


0.274 


Pv (m) 


0.596 


Dispersion* r]x (m) 


0.0 


Betatron tune 


13.368 


Uy 


10.398 


Dynamic aperture* Ux (mm) 


±1.75 


ay (mm) 


±3.50 


Acceptance ex (mm - mr) 


11 


ey (mm - mr) 


21 


Momentum acceptance (%) 


±1.2 


Number of longitudinal buckets 


203 



*At the center of straight sections 



The beta functions and the horizontal dispersion are 
shown in Fig. ^p. Small beta functions in the bends pro- 
duce a smaller beam profile, allowing the bend elements 
to be stronger and the beam to occupy the most linear 
region of the elements. The straight sections are designed 
to be free of horizontal dispersion to avoid synchro-beta 
coupling at the bunchers. 

If uncorrected, the vertical closed orbit displacement 
caused by gravity is 2.6 mm and is large enough to cause 
loss of the circulating beam. The orbit is corrected by 
displacing Q2 by 0.24 mm downward to produce upward 
kicks. The resulting vertical orbit distortion shrinks to 
0.26 mm as shown in Fig. [Sjand is no longer a problem. 

With this orbit correction, the dynamic aperture for 
400 turns, at the center of a straight section, is about 
2 mm by 3 mm half-width as shown in Fig. 2| This 
dynamic aperture corresponds to acceptances of 11 mm- 
mr horizontal and 21 mm-mr vertical, as listed in Table 
nil The resulting beam size is shown in Fig. Oi. The 
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FIG. 3: Beam half-widths (a) and the beta functions and 
horizontal dispersion (b) in the storage ring. Beta is the dis- 
tance in which the transverse (betatron) oscillation advances 
in phase by one radian. A schematic of the lattice is shown 
for location reference. 
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FIG. 4: Starting coordinates in the center of the straight 
section for the molecules that survive 400 turns. This defines 
the dynamic aperture. 



momentum acceptance, calculated by the multi-particle 
tracking simulation, is ±1.2% which is equivalent to an 
energy acceptance of ± 45 mK. 
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FIG. 5: Corrected vertical closed orbit displacement of the 
beam in the storage ring 



IV. DECELERATED BEAM 



TABLE III: Parameters of the decelerator for injected beam 



T ^ TTl pi" PV 
J. didiiiCLCi. 




\/p1 /"ir*! t^r ct'i' ci^iiT"r*o itti /c i 


■m n 

O-LyJ 


x/plr^r*! t^r a"f pvii~ im /o i 
VdUUlLV cLlf "AIL 1 111/ b ) 






ZIZZ 


Length, of bunch at exit (mm) 


1 n 

lU 


Tr.TYii"H" ci Ti /^PG Eli" P"V"ii" V CI n ri t7 rrnrn mvl 
J_Jllii L LdiiCcb cLL cAl L ^ jV cLiiU V I iliili ilil 1 


ou 


Electrode gap (mm) 


7 


Decelerating field at entrance (MV/m) 


9 


Decelerating field at exit (MV/m) 


4.5 


Length of last decel. electrode (mm) 


24 


Length of decelerator (m) 


19.6 


Number of decel. electrodes 


139 



A. Decelerator 

To reduce the velocity from the 310 m/s at the source 
to 30 m/s requires many stages of deceleration by pulsed 
electric fields in a long linear array. At each of the 139 
decelerating stages, a bunch of molecules enters a set of 
parallel electrodes when the field is zero; the field pulses 
on and the molecules lose kinetic energy equal to |Vl^(i?)| 
as they exit the electrodes. 

Our decelerator desim^differs in almost every way from 
previous designs 0, A decrease in the strength 

of the electric field while the bunch exits the electrodes 
provides longitudinal restoring action that prevents the 
bunch lengthening due to velocity spread |l5j| . The 
lengths of successive electrodes decreases as the velocity 
and spacing of the bunches decreases. 

Interspersed between the pulsed parallel electrodes 
are alternating-gradient lenses to confine the molecules 
transversely. Their overall focusing action must be 
stronger in the plane of the electric fields to counter the 
dcfocusing from fringe fields. The major parameters of 
the decelerator are summarized in Table Hill Details of 
decelerator design will be published later. 

B. Injector 

To inject the beam, we need a bend electrode that can 
pulse on or off in the time between buckets in the ring. 
This allows us to store multiple (up to 203) bunches in 
the ring. The deflecting electrode (Fig^ is part of a 
transport line that transforms the pulse leaving the de- 
celerator to match the orientation of the transverse ac- 
ceptances of the ring at the point of entry onto the closed 
orbit of the ring. The deflecting electrode is actually an 
array of bend electrodes with radius 0.6 m, similar to a 
bend section in the ring. A horizontal phase advance of 
2TT in this bend, avoids a net dispersal of molecules that 
are within the ±2 % velocity spread. 

In passage along the line, the velocity spread of ±2% 



lengthens the bunch and a debuncher at the point of in- 
jection (FigP) brings 90% of the bunch within the ±1.2% 
longitudinal momentum acceptance of the ring. 



C. Source and Intensity 

We calculate the intensity based upon a pulsed jet 
source with 1% methyl fluoride seeded in xenon carrier 
gas, using the equations in Miller jT^ and verified against 
seeded xenon jet source performance reported in the liter- 
ature 0,|3,[i3|- Xenon's high mass (133) produces much 
slower beams (310 m/s from a room-temperature reser- 
voir) than do light carrier gases, resulting in a shorter 
(19.6 m) decelerator. 

The bunch intensity is determined by the source fiow 
rate, the J — state population, the velocity distribu- 
tion and the acceptances. A source orifice of 1 mm di- 
ameter and reservoir pressure of 6.56 x 10^ Pa (500 Torr) 
will produce an intense cold beam with a peak intensity 
of 3 X 10^^ molecules sr~^ s~^, a longitudinal velocity 
spread of 7.2 m/s FWHM, and less than 1% clusters. 
We estimate the methyl fluoride J = rotational state 
fraction to be 30%. In an apparatus with a finite pump- 
ing speed, this peak intensity is possibly only by using a 
pulsed jet source operating with a small duty cycle. The 
short widely-spaced beam pulses entering the decelerator 
(which become more closely spaced after deceleration) re- 
quire a duty cycle of less than one percent for a 100 Hz 
pulse rate. This would allow all 203 buckets in the ring 
to be filled in 6.4 turns. 

The transverse and the longitudinal emittances (units 
of m^ s^^) of a bunch of molecules are unchanged in pass- 
ing through the deceleration process ^ISJ from the source 
to their injection in the storage ring. Therefore the frac- 
tion of molecules from the source that enters the ring is 
the product of the ratios of ring acceptances to source 
emittances. In the transverse directions, the beam from 
the source has ±0.5 mm spatial extents and ±1000 mr an- 
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gular divergences; then the horizontal and vertical accep- 
tances of the storage ring (Table ^ of 11 mm-mr and 21 
mm-mr respectively, result in 8.66 x 10"^ of the molecules 
being transversely accepted. Longitudinally, one second 
of beam from the source is 310 m long and has a velocity 
spread of ±3.6 m/s. The storage ring will accept ±0.6 
m/s in a 10-mm long bunch, which is 5.4 x 10~® of the 
source longitudinal emittance. 

Combining these nunbers and acounting for the 90% 
acceptance of the storage ring from the injector yields an 
intensity of 3.8 x 10* molecules/bunch. Bunches could 
be injected into the storage ring singularly or in large 
numbers. With a maximum of 203 stored bunches there 
would be nearly 10^^ molecules circulating in the storage 
ring and a flux of 2.5 x 10^^ molecules/s. Each bunch 



would have a density of about 3 xlO^ molecules/cm^ in 
the long straight sections, and higher in the bends. 



V. ACKNOWLEDGMENTS 

The authors acknowledge and thank Richard Gough 
and David Robin for their enthusiastic encouragement, 
and Swapan Chattopadhyay and Ying Wu for early con- 
tributions to the storage ring work. Work supported by 
the Director, Office of Science; Office of Basic Energy Sci- 
ences, and Office of High Energy and Nuclear Physics, 
U.S. Department of Energy, under Contract No. DE- 
AC03-76SF00098. 



[1] J. D. Weinstein, R. deCarvalho, T. Guillet, B. Friedrich 

and J.M. Doyle, Nature 395, 148 (1998). 
[2] H.L. Bethlem, G. Berden, F.M.H. Crompvoets, R.T. 

Jongma, A. J. A. van Roij, and G. Meijer, Nature 406, 

491 (2000). 

[3] F.M.H. Crompvoets, H.L. Bethlem, R.T. Jongma, and 

G. Meijer, Nature 411, 174 (2001). 
[4] F.M.H. Crompvoets, H.L. Bethlem, J. Kiipper, A.J.A. 

van Roij and G. Meijer, Phys. Rev. A69, 063406 (2001). 
[5] H. Nishimura, G. Lambertson, J. G. Kalnins, and H. 

Gould, Rev. Sci. Instr. 43, 3271 (2003). 
[6] A. Volpi and J.L. Bohn, Phys. Rev. A65, 052712 (2002). 
[7] J. L. Bohn, Phys. Rev. A63, 052714 (2001). 
[8] M. Kajita, T. Suzuki, H. Odashima, Y. Moriwaki, and 

M. Tachikawa, Jpn. J. Appl. Phys. 40, L1260 (2001). 
[9] M. Kajita, Eur. Phys. J. D20, 55 (2002). 
[10] A.V Avdeenkov and J.L. Bohn, Phys. Rev. A66, 0052718 

(2002). 



[11] J.G. Kalnins, G. Lambertson, and H. Gould, Rev. Sci. 

Instr. 73, 2557 (2002). 
[12] H. Nishimura, G. Lambertson, J.G. Kalnins, and H. 

Gould, IEEE Proc. of PAC 2003, 1837 (2002). 
[13] H.L. Bethlem, A.J.A. van Roij, R.T. Jongma and G. Mei- 
jer, Phys. Rev. Lett. 88, 133003 (2002). 
[14] M.R. Tarbutt, H.L. Bethlem, J.J. Hudson, V.L. Ryabov, 

V.A. Ryzhov, B.E. Sauer, G. Meijer, and E.A. Hinds, 

Phys. Rev. Lett. 92, 173002 (2004). 
[15] J. A. Maddi, T.P. Dinneen, and H. Gould, Phys. Rev. 

A60, 3882 (1999). 
[16] D.R. Miller, in Atomic and Molecular Beam Methods, 

edited by G. Scoles (Oxford University Press, New York, 

1988), Vol. 1, p. 14. 
[17] M. Gupta and D. Herschbach, J. Phys. Chem. 103, 10670 

(1999). 

[18] G. Lambertson, private communication, 2004. 



